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The directed transfer function (DTP) introduced by Kaminski and Blinowskal (Il99lh is a well- 
known frequ ency-domain based me asure for the interrelationships in multivariate time series. In 
the paper by iKamiriski et alj ( 2001 ). the authors claim a relationship between the DTP and the 
concept of Granger causality. Here, Granger causality from one channel Xi to another channel Xj 
is defined in terms of a bivariate VAR model 

Xi{t) = {2 Ai{u)Xi{t - u) + E A,j{u)Xj{t -u) + ei{t) 

u=l u=l 

Xj{t) = E Aji{u)Xi{t - n) + f: Ajj{u)Xj{t - u) + ej{t), 

u=l u=l 



, p. We note 



and Xi is said to Granger cause Xj if Aji{u) is nonzero for some u = 1, 
that this bivariate notion of Granger causali ty has been widely used (e.g.. Iplorens and Mouchart 
1985L iGoebel et al.ll2003L faesse et all 2003h . but for multivariate syster ns a more general notion 
of Granger causality in terms of multivariate VAR mod els exists (e.g.. ISims 1 19801 lHsiaol[l982L 
Toda and Ph ilipps"l993'.'Havo"l999'.'Eichleill2007Ll2005h. which is more in hne with the original 
definition by Granger (1969, 1980, 1988). As illustrated in Eichler (2006), the DTP is neither 
a measure for this multivariate Granger causality nor actually for the bi variate notion describe d 
above. In the following, we show in detail the problems with the proof of iKaminski et all (|200lh . 

Por the proof of a relation between bivariate Granger causality and DTP, the authors derive the 
bivariate autoregressive representation of two components of a multivariate VAR(p) process (cf 
eqs (12) to (14)). We note that the autoregressive representation of a weakly stationary process is 
defined in terms of linear projections, which implies that the error process e{t) = {ei{t), ej{t))' is 
white noise, that is, the errors at different time points are uncorrelated. In the frequency domain, 
this implies that the spectral matrix of the error process is constant and equal to S/27r, where 
S; = var(e(t)). 

In the paper, the authors derive the bivariate autoregressive representation (setting i = 1 and 
J = 2) expressed in the frequency domain 



[An(A) - Ai2(A)A22(A)-^A2i(A)] Tj^l^j 



E[{X) 



(1) 



(cf eqn (14)) with error process 



p3(A)^ 
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The spectral matrix fe^e'^W of the error e'{t) process is given by 



5^13 • • • ^Ip 
S23 • • • 



(A22(A)')-^Ai2(A)' 
/II33 ... Ssp^ 



+ Ai2(A)A22(A) 



-1 



(A22(A)')-'Ai2(A)' 



\^p3 • • • ^pp y 

Due to the frequency dependency of Aii(A), Ai2(A), and A22(A), this expression in general will 
not be constant over frequency and, thus, cannot be the spectral matrix of a white noise process. 
Consequently, the process e'{t) = (ei(t), e2(t))' defined by E'{X) = {E[{X) , E'2{X)) in general 
is not a white noise process and ^ is not the desired bivariate autoregressive representation. 

That e'{t) = (ei(t), e2(t))' indeed is not generally a white noise process can be shown by a 
simple example. Consider a simple trivariate VAR(l) model 

Xi{t) = aX3{t-2)+ei{t), 

X2{t) = ^X3it-l)+e2it), 

where e{t) = {£i{t),e2{t),e3{t)) is a white noise process with mean zero and variance equal to 
the identity matrix. On the one hand, we have 

/I 
A(A) =01 
\0 

and simple manipulations show that 

H(A) = A(Ari = 

which implies that the DTF from channel 2 to channel 1 is zero. 

On the other hand, the bivariate autoregressive representation is given by the best predictor of 
X{t) = (Xi(i),X2(t)) based onX{t- l),X{t - 2), . . .. It can be shown that it is given by 

X2it)=e2{t), 
where £2{t) = £2{'t) + /3^3(i — 1) and 

= 61 (t) - jf^ e2{t-l) + ^ es{t - 2). 

Note that = {£i{t),£2{t)) is indeed a white noise process satisfying 

E(e(t)e(s)') = 

for alH 7^ s. In particular, we have 

cov(.-i(t - l),e2{t)) = + = 0- 

It follows that X2 bivariately Granger ca uses X-j despite the fact that the DTF is zero. Thus the 
example contradicts the result by .Kaminski et al.i 




We note that the error proce ss e in the above bivariate representation differs from the error 
process e' proposed by .Kamiriski et all , which is of the form (written in the time domain) 

e'lit) = Slit) + ae^{t - I) 

e'2{t)=£2{t)+^£^{t-2)' 

Obviously we have 

E(e;(t-l)4(t)) = a/3^0, 
that is, the process e'{t) = {e'^{t) , e'2{t)) is not a white noise process as required by the autore- 
gressive representation used in the definition of Granger-causality. As a consequence, the temporal 
dependence structure that is still hidden in the dependencies of e' is neglected when computing 
Granger-causality based on the bivariate representation ([T]). 
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